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Abstract 

We investigate the double Dalitz decays B s — > l + l~l' + l'~ on the basis of 
the effective Hamiltonian for the transition bs — * l + l~, and universal form 
factors suggested by QCD. The correlated mass spectrum of the two lepton 
pairs in the decay B s — ► e + e~ is derived in an efficient way, using a 

QED result for meson decays mediated by two virtual photons: B s — > 7*7* — ► 
e + e~ . A comment is made on the correlation between the planes of the 
two lepton pairs. The conversion ratios piwi' = ^(b ^77) — are es timated 
to be p eeee = 3 x 10~ 4 ,p eew = 9 x 1CT 5 and = 3 x 1CT 5 , and are 

enhanced relative to pure QED by 10 — 30%. 
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1 Introduction 



In a recent paper [T] we investigated the decay B s — > = e,/x), using the 

effective Hamiltonian for the transition bs — > , and obtained a prediction for the 
conversion ratio 

p " - r(5 s 77) (1) 

in terms of the Wilson coefficients CV, C9 and C\q. An essential ingredient of the 
calculation was the use of a universal form factor characterising the matrix elements 
(j\si<j^ u (l +7 5 )6|-B S ) and (7(57^(1 ± 75 )6|i? s ), as suggested by recent work j2] on 
QCD in the heavy quark limit (m& 3> Aqcd)- It was found that the ratio pu 
was significantly higher than one would expect from a QED calculation of Dalitz 
pair production B s —>■ 7*7 — > Z + Z~7, the difference reflecting the presence of the 
short- distance coefficients C$,Cio, as well as the universal behaviour of the 

QCD- motivated form factor. The purpose of the present paper is to apply the same 
considerations to the "double Dalitz decay" B s — > , to determine whether 

there is similar enhancement of the double conversion ratio 

T(B 8 -+i+rv+i'-) 

Plll ' V - T(B S -+ 77) ' { ) 

compared to what one would obtain from the QED process B s — > 7*7* — > . 
We examine also the correlation in the invariant mass of the two lepton pairs, and 
the nature of the angular correlation between the and Z /+ Z'~ planes, which is a 
crucial test of the B s — > 77 vertex. 

2 Matrix Element and Invariant Mass Spectrum 

We begin with the effective Hamiltonian for bs — > l + l~ jS] 

n cS = ^v tt v: s \cf{s^p L b)h»i 

+ C 10 {s lfl P L b)l W l (3) 

c 7 _ 

- 2—sia lJlu q v (m b P R + m s P L )bl-i^l 
q 2 



where Pl.r — (1 ^75)72 and q is the sum of the Z + and I momenta. Ignoring small 
q 2 — dependent corrections in Cg S , the values of the Wilson coefficients are 

C 7 = -0.315, C 9 = 4.334, C 10 = -4.624 . (4) 
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Then, as shown in jl], the matrix element for B s — > l + l 7 has the form 



V2n m ts M Bs 
■ [e^pa^q'k^Aj^l + A 2 h^ 5 l) (5) 
+ i{e*{k ■ q) - (e* ■ qJk^Bj^l + B 2 lj^ 5 l)] 



where 



M Bs 



A 2 — Cio/v 



B\ = Cgf'A + 2C7 ^-ffp , 

qi 

B 2 = C'io/a • 

The form factors fy, /a, fr, fxi defined in Ref. pQ, will be taken to have the universal 
form 

f r = f A = fT = f I= l^l. + eK ^), (7) 

predicted in the heavy quark approximation (ra& S> Ag C D, m b ^> m s) in QCD j2]. 
Here, A s = m Bs - m b « 0.5GeV,x 7 = 2E 1 /M Bs = 1 - ? 2 /M| s , and / Bs « 200MeV 
is the 5 S decay constant. The essential feature for our purpose will be the universal 
1/x-y behaviour, the absolute normalization dropping out in the calculation of the 
conversion ratio. (Corrections to universality are discussed in Ref. jS]). 

To obtain the matrix element for B s — ► we treat the second lepton pair 

as a Dalitz pair associated with internal conversion of the photon in B s — > 
/ + /~7. From this point on, we will specialise to the final state e + e~/i + /i~, consisting 
of two different lepton pairs. This avoids the complications due to the exchange 
diagram that occurs in dealing with two identical pairs. The matrix element then 
has the structure 

M(B S - e+e-^V) ~ ^{a + {q 2 )L%{ qil q 2 ) + a-(q 2 )L»( qi , q 2 ))L» m (h, k 2 ) 

K (8) 

• [e^ P aq p k a + i(g^ u k ■ q - k^q u )] 

where k and q are the four-momenta of the two lepton pairs, k 2 and q 2 being the 
corresponding invariant masses. The currents L± and L em are given by 

LUQUQ2) = u(qi)Y(l ±l 5 )v(q 2 ) , 
L» m (k 1 ,k 2 )=u(k 1 ) 1 > 1 v{k 2 ). [ j 
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where k% + k 2 = k,q% + q 2 = q. The coefficients a±(q 2 ) are related to those in Eq. 
® by 

a ± (q 2 ) = AM 2 ) ± A 2 (q 2 ) , (10) 

where we have used the fact that for universal form factors, B\ 2 = A\ t2 . 

At this stage, it is expedient to compare the matrix element (JHJ) with the ma- 
trix element for double Dalitz pair production in QED. We will make use of the 
recent analysis of Barker et al. [0], who have studied the reaction Meson — > 7*7* — > 
, using a vertex for Meson — > 77 that is a general superposition of scalar 
and pseudoscalar forms, the matrix element being 

e e 

k 2 g 2 — «... - ^ 

■ [fp + €s($iu>k ■ q - k^q u )} . 



M Barker = const. ■ — — L% m (qx, q 2 )L" m (ki, k 2 ) 



The coefficients £p and £s are normalized so that |£p| 2 + l^l 2 = 1. (In Ref. [0j they 
are denoted by £p = cos £,£5 = sin£e i<5 .) 

From this matrix element, Barker et al. have derived the correlated invariant 
mass spectra for the decay into e + e" / u + 'fj,~ (ignoring form factors at the M7*7* 
vertex) 



1 



77 



d 2 T 

dxi 2 dx34 



- Barker 



9tt 



A|2A:!4A (3-A 2 2 )(3-A 2 4 ) 



l£p| 2 A 2 +fc| 2 (A 2 + ^) 

The variables entering the above formula are defined as follows: 

X12 = (qi + g 2 ) 2 /M 2 = q 2 /M 2 , 

X34 = (h + k 2 f/M 2 = k 2 /M 2 , 
2 

mi 



(12) 



^ 2 1 



Xi = x 2 



M 2 x 



12 



m 



X3 
z 
A 12 

A34 

w 2 



X4 
1 - 



M 2 x u ' 
X12 — X34 . 



(13) 



a/(1 - xi - x 2 ) 2 - 4xix 2 
\/(l — X3 — X4) 2 — 4x3X4 

4^12^34 , 



A = 



W7 Z 



4 



Here mi and denote the masses of the electron and muon, and M the mass of the 
decaying meson. The phase space in the variables x\ 2 and x 34 is defined by x® 4 < 



34 

x 34 < (1 - y/xu) 2 , £? 2 < x i2 < (1 - y/xu) 2 , where x° l2 = Amf/M 2 , x® u = 4m 2 i /M 2 . 

We can now adapt the QED result (|T2j) to the process B s — > e + e~fi + fi~, by 
comparing the matrix element (JTTj) with that in Eq. (jHJ). The essential observation 
is that in the approximation of neglecting lepton masses, the vector and axial vector 
parts of the chiral currents contribute equally and independently to the invariant 
mass spectrum. In addition, the matrix element for B s decay corresponds to the 
QED matrix element considered by Barker et al., if we put £p = l/v2, £s = i/v2. 
This allows us to obtain the invariant mass spectrum for the double Dalitz decay 
B s — > e + e~ in electroweak theory: 



. 77 



dT 



dx^dxu 



EW 



(V9 + — ) 2 + v!o 
X\2 



+ 



™,2 ,y,2 

x 12 x 34 

<->-» 2 1 ,7-* 2 
x 12 ^ x 34 



|^(a;i2,£34) 



XZA 

dT 



T~~ yC?Xi2'^2'34 



(14) 



QED 



where 



L 77 



dxi 2 dx 34: 



QED 



a 2 A12A34A 

97T 2 W 2 



(3-A 2 2 )(3-A 2 4 )(2A 2 + V). (15) 



Here we have used the abbreviation 779 = Cg/(2Cj) and 7710 = C\q/(2Ct), introduced 
in Ref. pQ. The electroweak formula (|T4*|) reduces to the QED result in the limit 
V9 = Vw = 0, F(x 12 , x 3i ) = 1. 

The form factor F(xi 2 , X34) is chosen to have the universal form 



^12, £34) 



xu 



X34) 



(16) 



(a possible normalization factor drops out in the calculation of the conversion ratio). 
This is a plausible (but not unique) generalization of the universal QCD form factor 
1/(1 — £12) that occurs in the single Dalitz pair process B s — > e + e~ r y. 

In Fig. ^we plot the correlated invariant mass spectrum for B s — > e + e~fi + fi~ 
in electroweak theory. The ratio of the electroweak and QED spectra is shown in 
FigEl and indicates the effects associated with the coefficients 779 and r] W , and the 
form factor F(xi2, 234). One notes a slight depression in the region x\2 = — 2 §f or 
X34 = — 2 §^-, connected with the vanishing of the term (C 9 + f^j) 2 or (C 9 + f^ 1 ) 2 - 
There is also a general enhancement for increasing values of Xi 2 , 234, because of the 



5 



form factor If the form factor F{xi2-,x 34 ) is set equal to one, the ratio of the 

electroweak and QED spectra has the structure plotted in Fig. EU illustrating the 
effects which depend specifically on the electroweak parameters rjg, r] W . 

The absolute value of the conversion ratio p ee ^ is obtained by integrating 
(j^— dT / dxi2<ix34) ew over the range of xu and 234. In the QED case, this ratio is 
conveniently expressed in terms of the integrals J1...6 introduced in Ref. jHj: 

I i = t;[ [ dx 12 dx 34 ^^Y— > 



3 J J w 



h = H I dx l2 dxJ lAXz2 



2 



3 J J w 
h = - / dx 12 dx u - 



3 / / w 



2 



I4 — I I dx^dx^ — (3 — \\ 2 — A34) 




(17) 



h = J J dxi 2 dx 3 ^ 12 ^ Z (3 - \\ 2 - A34) , 
h = g J J rfxi 2 rfx34Ai2A 3 4A(3 - Ai 2 )(3 - Ag 4 ) . 

These integrals are listed in Table Q (where, for completeness, we have also given 
the values for the final states eeee and /jfififi). These integrals allow us to calculate 
the QED double conversion ratio 

a 2 



(Peew)QED = ^2~( J 1 + J 2 + 2(J 4 + I 5 + J 6 )) ^ 

= 7.6 X 1(T 5 . 

The corresponding result for electroweak theory, based on the differential decay rate 
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(|14|) . can be expressed in terms of the integrals 
h 



h 
h 
h 
h 




A 3 A 3 A 3 
dx 12 dx 3 4 12 - 3 4 — G(xi2, x 34 ) 





w- 

dx^dx^ - G(xi2, X34) 

W 

dXi2dx%4 G(Xi2, X34) 





dxi 2 dx 34: 
dx^dx^ 



itr 

A12A34A 3 



(19) 



A12A34AZ 



(0- 



(3 — A^ 2 — A3 4 )G(xi2, x 34 ) , 
2 

-(3 — A 12 — A 34 )G(xi2, x 34: ) 




c?xi 2 (ix34Ai2A34A(3 — A 42 )(3 — A3 4 )G(xi 2 , £34) 



The factor G(a;i2, X34) in the integrand of Eq. (fT9| contains the effects of the elec- 
troweak coefficients 779, 7710 and the universal form factor -F(x 12 , £34): 



G(Xi2, X34) 



V9 



X12 



Vw 



+ 



V9 



£34 



Vw 



„2 ^2 



(20) 



a -12 X 34 1 T7if„ ™ M2 



,71 2 I ,-y»2 

x 12 "r ^34 



\F(x 12 ,X U ) 



The integrals Ii,...,Iq are given in Table El The electroweak conversion ratio, 
analogous to the QED result (fT8|) . is given by 



(p. 



a 



eeiifij EW 



6ir 2 
9.1 x 10~ 5 



h + h + 2(/ 4 + h + h 



(21) 



In comparison to the QED result ()18|) . the double conversion ratio for B — > ee/i/i in 
electroweak theory is enhanced by ~ 20%. 

A calculation of the spectra for the channels eeee and fifififi is complicated by 
interference between the exchange and direct amplitudes. The conversion ratio for 
these channels takes the form 



P = Pi + P2 + P\2 



(22) 
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where p\ and p 2 denote the "direct" and "exchange" contribution, and p i2 an in- 
terference term. Numerical calculations of the decays 7r° — > e + e~e + e~ and Kl — > 
e + e~e + e~ suggest that pu is small and p\ ~ p 2 . Thus a rough estimate of the double 
conversion ratio can be obtained using the formula (|2ip. with an extra factor (i) • 2 
where (|) is the statistical factor for two identical fermion pairs, and 2 comes from 
adding direct and exchange contributions. This yields, using the numbers in 
Table H 

{Peeee)EW « 2.9 X 10~ 4 , 

(/WW « 2 -8 x iO^ 5 • 
For comparison, the QED results, using Tabled are 

(p ee - ee -)Qi?D ~ 2.7 x 10~ 4 , 
(p»W~ 2 -2x 10 -5_ 



Thus the enhancement in the case of eeee is ~ 10% and that in pppp about 30%. 
Combining (21) and (23), the ratio of the channels eeee, eepp and pppp is approx- 
imately 

eeee : eepp : pppp 
= 3 : 1 : 0.3 ^ ' 

To obtain the absolute branching ratios, we note that the decay rate of B s — > 77, 
derived from the effective Hamiltonian (jSJ), involves the Wilson coefficient CV and 
the universal form factor /r(^ 7 = 1) (see Eq. (J7|)). Using nominal values for fs s 
and A s , and evaluating CV at the renormalization scale p = mj, Ref. j7j finds 
Br(5 s — > 77) = 1.23 x 10 -6 . Using this as a reference value, we obtain: 

Br(B s -> eeee) = 3.6 x 10~ 10 , 

Br(B s -> ee/ip) = 1.1 x 10~ 10 , (26) 
Br(5 s -> ppp/i) = 3.5 x 10~ n . 

3 Correlation of e + e~ and /j, + fi~ planes in B s —> 

One of the distinctive features of the electroweak B s — ► 77 matrix element is that 
the coefficients £s and £p (normalized to \^s\ 2 + |£p| 2 = 1) are given by £5 = 
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A= and = The equality |£s| 2 = |£p| 2 leads to the simplification that the 



factor |6>| 2 A 2 + |6| 2 (A 



2 + %w 2 ) appearing in the spectrum (fEZJl could be written as 



|[2A 2 + §u> 2 ] in going over to the electroweak case (Eq.(15)). A further interesting 
consequence is the distribution of the angle between the e + e~ and planes in 

B s — > ee/i/2. Generalising the QED result given in Ref. [U] to the electroweak case, 
the correlation in is given by 



i dr\ ee ^ 

77 d4> J ew 



a 



67T 3 - 



h sin 2 + J 2 cos 2 + (J 4 + h + h) 



(27) 



The fact that I 2 is so close to I\ means that the spectrum dT/dip is essentially 
independent of 0. Furthermore, the fact that a>rg(£s/£p) = 7r/2 refelects itself in 
the absence of a term proportional to sin cos 0, the presence of which would lead 
to an asymmetry between events with sin cos > and < 0. 

It may be remarked that there are corrections to the B s — > 77 matrix element 
(associated, for example, with the elementary process bs — > cc — >• 77) which cause 
the superposition of scalar and pseudoscalar terms to deviate slightly from the ratio 
£,s/£,p — i [3 IE]- From the work of Bosch and Buchalla [Zj, we find 



' 2C 1 +NC 2 \ B 

1 - o n 9 ^ c 

6 C7 VTiB 



(28) 



-2 + (-2 In 2 z + 2tt 2 



and z c = m 2 c /ml ~ 0.1, G\ 



where 

4iri\nz)z + 0(z 2 ) , (29) 

1.1, Gi = — 0.24, iV = 3. There is thus a small 
correction to the equality |£p| = More interestingly, the phase 5 = arg(£p/£s) 
is not exactly 90°, implying that a term of the form ^3 sin cos cos 5 could appear 
in dY/d<p. These corrections are, however, too small, to have a measurable impact 
on the spectrum and branching ratio of the decay B s — > ee\x\x calculated above. 



4 Conclusions 

We have calculated the spectrum and rate of the double Dalitz decay B s — > e + e~ , 
using the effective Hamiltonian for the flavour- changing neutral current reaction 
bs — > , and form-factors motivated by the heavy quark limit of QCD. A method 
is given for obtaining the correlated mass spectrum dY '/ ' dxyidx^ from the known 
results for the QED process B s — > 7*7* — > e + e~ / u + . The conversion ratios p^yy = 
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r(B s — > l + l~l' + l'~) /Y{B a — > 77) show an enhancement over the QED result, ranging 
from 10% for the channel e + e~e + e~ to 30% for the channel n + n~ n + n~ . Our best es- 
timate of the branching ratios, using the QCD estimate By(B s — ► 77) = 1.23 x 10~ 6 
given in fT], is Bi(B s -> eeee) = 3.6 x lO" 10 , Br(5 s -> ee/i/i) = 1.1 x lO" 10 , Br( J B s . -> 
HjxHp) = 3.5 x 10 -11 . These branching ratios may have a chance of being observed 
at future hadron machines producing up to 10 12 B s mesons. 
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B s — > ee/jpL 


B s — > eeee 


B s — > /i/i/X/X 


h 


7.754 


32.501 


1.772 


h 


7.806 


32.556 


1.821 


h 


15.556 


65.053 


3.589 


h 


17.558 


58.416 


5.115 


h 


17.641 


58.499 


5.199 


h 


0.0548 


0.0556 


0.0540 



Table 1: Numerical values of the integrals I\, . . . , J 6 for 5 S -> ///'/' in QED 





B s — > ee/i/i 


B s — > eeee 


B s — > /i/i/i/x 


h 


9.336 


35.491 


3.856 


h 


9.477 


35.643 


4.002 


h 


18.793 


71.114 


7.837 


h 


20.411 


63.457 


5.784 


h 


20.637 


63.685 


6.003 


h 


0.148 


0.152 


0.146 



Table 2: Numerical values of the integrals I±, . . . , J 6 for £> s — > in electroweak 
theory 
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Figure 1: Invariant mass distribution dY/dx^dx^i for B s — > e + e in elec- 

troweak theory 
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Figure 2: Ratio (dT / ' dx^dx^) ew I '(dT / ' dx^dx^QED showing influence of elec- 
troweak parameter C 9 , C w and the form factor F(xi 2 , x 3i ). 
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Figure 3: Ratio (dT / 'dx^dx^EW I '{dT / ' dx^dx^QED in the limit of a constant form 
factor (F(xi2,x^i) = 1), illustrating specific effect of electroweak parameters rj 9 = 
C 9 /(2C 7 ) and rj 10 = C 10 /(2C 7 ). 
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